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Abstract:  We  study  the  stability  with  probability  one  of  the 
stochastic  bilinear  system  dX  =  AX  ds  +  BX  dw,  where  A  and  B  are 
fixed  matrices  and  w  is  a  Brownian  motion.  Bounds  for  the 
Lyapunov  numbers  associated  with  this  equation  are  given. 

Bilinear  noise  models  are,  after  linear  ones,  the  second  simplest 
case  of  stochastic  systems;  they  may  arise  in  many  problems  in 
which  linear  noise  models  are  inappropriate  (many  examples  are 
given  in  [6]). 

The  aim  of  this  paper  is  to  give  a  condition  for  the  stability 
with  probability  one  of  the  d-dimensional  Ito  equation  which 
describes  the  behavior  of  such  a  system 

dYs  =  AYs  ds  +  BYsdws  (1) 

Yo  =  y 


where  A  and  B  are  two  given  dxd  matrices  and  w  is  a  scalar 
Brownian  motion  (see  also  the  more  general  equation  (12)  below). 
I.e.,  we  want  to  find  whether  or  not  Ys  tends  a.e.  to  zero  as  s  tends  to 
infinity.  Note  that  in  the  one  dimensional  case,  we  have  an  explicit 
solution 

Ys  =  Yo  exp{t(A-B2/2)  +  Bwt), 

and  the  stability  with  probability  one  is  guaranteed  iff  2A-B2<  0.  Note 
also  that 

E[  (Ys)2  ]  =  (Yo)2  exp2t(A+B2/2) 
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and  the  stability  is  guaranteed  iff  2A+B2  <  0.  This  is  a  much 
stronger  condition  whose  generalisation  in  d  dimensions  is  easy, 
thanks  to  the  following  equations: 

dYs®Ys  =  (I®A  +A®I  +  B®B)Ys®Ys  ds  +  (I®B  +  B®I)Ys®Ysdws 

dE[Ys®Ys  ]=  (I®A  +A®I  +  B®B)  E[Ys®Ys]  ds 

(for  the  definition  and  the  basic  properties  of  ®,  see  [l]).The  L2 
stability  is  then  governed  by  XmaxCI  ®A  +A®I  +  B®B)  (in  this  paper, 
^max(M)  denotes  the  largest  real  part  of  eigenvalues  of  the  matrix  M, 

and  Amin(M)  =  -Xmax(-M)  ). 

We  will  give  an  upper  bound  for  the  largest  Lyapunov  number 
of  (1)  (the  smallest  y  satisfying  the  bound  of  theoreml).  Note  that 
the  only  existing  expressions  for  the  largest  Lyapunov  number  of  (1) 
or  (12)  with  general  matrices  are  asymptotic  expansions  (B=8Bo,  e 
tends  to  0,  A  fixed)  in  dimension  2  ([5]).  The  following  criterion  will 
be  proved  : 

Theoreml.  Setting 

y  =  “  Vax(I  ®(A-B2)  +(A-B2)®I  +  B®B), 

then,  for  any  value  of  Yq, 

limsup  ^log(IIYtll)  <  y  a.e. 

t-^  *■ 


Considering  the  matrix  equation 

dPs  =  APs  ds  +  BPsdws 

(2) 

Po  =  I, 

Ys  can  be  expressed  as 


Ys  =  PsYo 
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and  the  following  basic  theorem  about  Lyapunov  numbers  (the  Xi's 
below)  is  true: 

Theorem2([2]):  There  exist  d  real  numbers  ....S  such  that 


lim  {PtT(co)Pt(co) }  i/2t  =  OT(o))AO((o)  a.e. 

t — 

where  O(co)  is  an  orthogonal  matrix  and  A  is  the  diagonal  matrix 
satisfying 


Ai,i  =  Qxp(Xi). 

Furthermore,  if 

Ei  =  span{0'r(a))ei,  OT(a))ei+i, . OT(o))ed  } 

(  (ei)i=i,d  is  the  canonical  basis  of  then,  a.e.,  we  have 
VueEi(co)\Ei+i(o)),  lim  -  LogllPt(co)ull  =  X,i. 

t-^oo  * 

All  the  are  equal  if  and  only  if  there  exist  a  matrix  M  such  that 

1,1, 

the  matrices  M(A  -  ^  trace(A))M'i  and  M(B  -  ^  trace(B))M'i  are 

skew-symmetric. 

Note  that  the  Lyapunov  numbers  of  Qs=Ps'^  are  (-^d."”»-^i)- 
To  prove  theorem  1,  instead  of  trying  to  find  upper  bounds  on  the 
Lyapunov  numbers  of  Ps,  we  will  get  lower  bounds  for  the  ones  of  Qs. 
Theorem2  shows  that  the  smallest  Lyapunov  number  is  generally 
attained  on  a  random  1 -dimensional  space;  this  fact  is  at  the  origin  of 
all  the  difficulty  of  stepl  below.  The  motivation  of  this  method  will 
appear  in  step2,  when  the  (non-quadratic)  term  (u'^Fu)2  in  (3)  will  be 
bounded  from  below  by  u^FTFu  (dealing  with  Pg,  we  would  have  to 
find  an  upper  bound  for  (uTBu)2), 

Proof  of  theorem  1 

c  denotes  a  constant  depending  only  on  d  and  is  allowed  to  change, 
finitely  often,  during  the  calculation. 

Note  that  Qs=Ps'^  satisfies  the  equation 
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dQs  =  EQs  ds  +  FQgdwg 
where 

E  =  -AT  +  B2T 
F  =  -BT. 


By  virtue  of  theorem2  (and  the  following  remark),  the  largest 
Lyapunov  number  of  Pg  is  a.e. 

-  lim  inf  “log(liQtxil). 

t^oo  X  ^ 

For  any  Xq,  Xt=QtXo  satisfies  the  equation 
dXg  =  EXg  ds  +  FXgdwg 

Using  Ito's  formula,  and  setting  Us=Xg/!iXgll,  we  obtain 

dxJXg  =  xJXs(  2usEus  ds  +  2usFus  dwg  +  u^FTFug  ds  ) 


dlog(IIXgll)  =  {  UgEug  +  ju]^FTFug  -  (ujFug)^  }ds  +  UgFug  dwg. 


This  equation,  due  to  Khas'minski  ([3]),  is  basic  for  the  study  of 
Lyapunov  numbers,  since  it  may  be  written 


t-i  {log(IIXtll) 


-  log(ilXoll)} 


t 


=t-i 


{usTEug+  jUgTFTFug- 

J 

0 


(ugTFus)2}ds) 


(3) 


t 

+  t‘i  JugTpUgdw 
0 


t 

stepl :  E[sup  I  I  UgTFugdwgl  ]=  o(t),  t  integer  >0. 

Xo  0 

The  sup  may  be  taken  a  over  a  dense  subset  of  Rd  (because  the 
stochastic  integral  is  a.e.  a  continuous  function  of  Xq)  so  that  the 
expectation  is  well  defined. 
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This  is  the  most  difficult  step  of  the  proof.  All  our  efforts  will 
be  focused  on  getting  rid  of  the  sup.  We  will  suppose  that  the  Xi  are 
not  all  equal  (if  they  are,  using  theorem2,  the  expression  above  is  wt 
times  a  constant  and  the  result  is  obvious). 

The  following  lemma  is  needed 

Lemma  Let  a  be  a  uniform  measure  on  the  the  unit  sphere  Sd-i  of 
X  any  vector  of  fsLi(c>),  P  an  invertible  matrix,  and  Q=¥-'^,the 
following  equalities  are  true: 

Cl  Jsign«x,y»  y  a(dy)  ci  =  (a(lyil))-i 

(4) 

J  f(y)a(dy)  =  det(P)  f(i||^)  llPylPd  (5(dy) 

(5) 

=  ci2jsign«x,y»sign«x,z»(zTpTFPy)  IIPyli-d-l||Pzll-d-l 

(6)  det(P)2  o(dy)a(dz) 

Proof  of  the  lemma: 

Observe  that  identity  (4)  has  only  to  be  proved  for  11x11=1,  and 
that,  because  of  the  rotational  invariance  of  a,  we  can  suppose 
x=(l,0,0....0)),  and  we  get  ci. 

The  left-hand  side  of  (5)  is  equal  to 

”  y,  r  Py 

""  cjf(jj^)  e-IIPxlldet(P)  dx 
=  c  det(P)  Jf(iil^)  e-rllPyll  rd-l  a(dy)  dr 

=  c  det(P)  llPyll-4  a(dy) 

the  constant  c,  independent  of  P,  is  identified  by  taking  P=I. 

In  order  to  prove  (6),  use  (4)  and  (5)  to  obtain 

ij^=  c  Jsign((x,QTy))  y  o(dy) 

=  det(P)  Jsign«x,y))Py  IIPyll-d-i  a(dy) 
which  implies  (6). 
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This  lemma  gives  us 
_  XqTQsTFQsXo  f 

UsiFus=  iio^X^lii  =  J  ^(Xo,y,z)'i'(Ps,y,z)a(dy)(7(dz) 
where 

d)(x,y,z)  =  sign«x,y»  sign«x,z» 

'P(P,y,z)  =  c  (zTpTpPy)  llPyli-d-iiiPzIl-d-i  det(P)2. 

And 


t  t 

JusFusTdws=  f  Jo(Xo,y,z)'P(Ps,y,z)a(dy)a(dz)  dws 

rv 


c  ^ 

0(Xo,y,z)  J^(Ps,y,z)dws  a(dy)a(dz) 


t 


sup  I  J  UsFus^^dwsI 
Xo  0 


r  t 


< 


J 


I  J'P(Ps,y,z)dwsl  a(dy)CT(dz) 
0 


t  r  t 

E^up  I  J  UsFus'^dwgl  ]  <  E[{  J'T2(Ps,y,z)ds  }  1/2  ]  a(dy)a(dz) 

Xo  0  Jo 


t  r  t 

Eiup  I  J  UsFus'Tdwgl  ]  <  E[  {  J'p2(Ps,y,z)ds  }l/2  a(dy)a(dz))  ] 

Xo  0  jo 

(7) 


we  will  prove  that  for  almost  all  (co,y,z), 

oo 

J'P2(Pg^y^Z)dS  <  oo 
0 

and  that  the  integrals 


(8) 
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|- (  J'p2(Ps,y,z)ds)l/2  (9) 

0 

are  uniformly  (on  t)  integrable  (over  (co,y,z)).  And  then  (7),  (8),  and 
(9)  will  imply  the  bound  of  stepl. 

To  prove  (8),  note  that 

l'P(Ps,y,z)l  =  c  (zTp^’pPsy)  IIPsyll-d-i|IPszll-d-i  det(Ps)2 

<c  liFII  IIPsyl|-d||Pszll-d  det(Ps)2 
<c  IIFII  exp{(-2?iid+2X?ti+  e(s))s} 


where  e(s)-^0  as  s-»<»  and  the  A,i  are  the  Lyapounov  numbers 
associated  with  P.  Since  we  are  in  a  case  where  the  Xi  are  not  all 
equal,  we  get  the  exponential  decrease  of  the  expression  above,  and 
(8)  is  true. 

To  prove  the  uniform  integrability  of  (9),  we  will  show  that 
r  t 


E[ 


a(  ^  ( J'P2(Ps,y,z)ds)i/2)  a(dy)a(dz)  ]  <  K 
0 


(10) 


where  K  is  a  constant  independent  of  t,  and 
a(x)  =  Ixl  logi/2(i+x2). 


Denoting 

t 

\l/(t,y,z)  =  Jl'P(Ps,y,z)|2ds, 
0 


equation  (10)  reduces  to 


-\|/(t,y,z)i/2  logl/2(l4-\i/(t,y,z)/t2)  o(dy)a(dz)  ] 


<  K. 


Stepl  will  be  finished  by  proving 


E[J  v(t,y,z)i/2  logi/2(l+\i/(t,y,z))  a(dy)a(dz)  ]  <  Kt. 

(11) 
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An  elementary  calculation  shows  that  the  function  xi/2iogi/2(i+x)  is 
concave;  this  implies 

(x+y)  l/2iogi/2('i+x+y)  <  xi/2iogi/2(l+x)  +  yi''2iogi/2(i+y)  for  any 

x,y  >0 


Defining  \|ri  by 

\l/i(t,y,z)  =  \ir(t+l,y,z)  -  >i^(l,y,z)  = 


t+1 

Jl'P(Ps,y,z)|2ds 

1 


we  obtain 

E[J  \|/(t+l,y,z)l/2  logl/2(l+\|/(t+l,y,z))  a(dy)a(dz)  ] 

<  E[J \|/(l,y,z)l/2  logi/2(l+\j/(l,y,z))  o(dy)a(dz)  ]  + 

E[J  \i/i(t,y,z)i/2  logi/2(l+\|/i(t,y,z))  o(dy)a(dz)  ]. 

We  will  now  get  the  index  1  out  of  the  last  formula.  The  basic  tool  is 
the  following  identity  resulting  from  Markov  property 


Ps+i  =  (0lPs)  Pi 


where  0i  is  the  shift  operator  on  the  trajectories  of  Brownian  motion. 
Using  this  identity  in 
t 

Vi(t,y,z)  =  JlT'(Ps+i,y,z)|2ds 
0 

and 

'P(P,y,z)  =  c  (zTpTpPy)  IIPyl|-d-l|lPzl|-d-l  det(P)2, 
we  get 

t 

Vi(t,y,z)  =  JlT^(eiPs,Piy,Piz)|2ds  det(Pi)2. 

0 


The  second  term  of  the  right-hand  side  may  be  rewritten  as 

E[J  \|r(t,Piy,Piz)i/2  logi/2(i4.xj/(t,Piy,Piz)det(Pi)4) 

det(Pi)2a(dy)a(dz)] 
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where  Pi  and  \i/(t,y,z)  are  independent(i.e.  constructed  from  two 
independent  Brownian  motions  and  the  expectation  is  taken  over  the 
product  space).  Using  lemma  1  we  obtain,  after  some  reductions  (Qi  = 


E[J  \}/(t,y,z)l/2  logl/2(i+x|/(t,y,z)det(Pi)4  liQiy||-2d  llQizih^d  ) 
a(dy)a(dz)] 

and  with  the  inequalities  (a,b>0) 

log(l+ab)  <  log(l+a)  +  log(l+b) 

(a+b)i/2  <  ai/2  +  bi/2 

we  get  the  upper  bound 

E[J  a(\|/(t,y,z)i/2)  a(dy)a(dz)  ]  + 

E[J  logi/2(i+det(Pi)4  llQiyl|-2d  ||Qizll-2d  )  a(dy)a(dz)] 

The  last  term  is  finite  since  log(IIQsyll)  satisfies  equation  (3)  and 
det(Ps)  satisfies 

dlog(det(Ps))  =  {trace(A)  -1/2  trace(B2)}  ds  +  trace(B)  dwg  ). 
Finally,  we  have 

E[J  a(\i/(t-i-l,y,z)l/2)  a(dy)a(dz)  ]  <  E[J  a(\i/(t,y,z)i/2)  cT(dy)a(dz)  ] 
-i-K 

where  K  is  a  constant.  This  ends  the  proof  of  (10),  and  the  first  step. 


step2:  end  of  the  proof. 

Using  stepl,  one  can  find  a  sequence  tn  such  that 


lim  —  sup  I  J  UsFus’^dwsI  ]  =  0  a.e. 

n— ^  Xq  q 


and  obtain,  with  equation  (3): 
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lim  inf  — log(IIXtnll)  >inf  u'rEu+ -uTfTFu-(uTFu)2. 
n->  oo  Xq  ‘•n  ueSd-i  ^ 


But 

(Fu  -  u(u"^Fu))T(Fu  -  u(u"^Fu))  >  0 
implies 

uTpTpu  >  (uTFu)2, 
which  gives  the  bound 


lim  inf  J-log(IIXtll)  >inf  uTEu  -  |uTfTFu 
t->oo  Xq  ^  ueSd-l  ^ 


>  X  +  inf  uT(E-^I)u  -  ^uTfTFu 
U€Sd-l  2 

for  any  real  X,. 

Actually,  we  have  still  the  choice  of  the  basis  of  Rd  ^ve  are  working 
in.  A  change  of  basis  with  a  matrix  R  changes  E  into  RER'i  and  F  into 
RFR‘i  and  then  for  any  invertible  matrix  R,  we  have 

lim  inf  J-log(IIXtll)  >  X  +  inf  uTR(E-XI)R-1u  -  ^uTR-TfTrTrfr- 

t->oo  Xq  ueSd-1  ^ 


lu 

>x  -  sup  vTS(-E+XI)v  +  ^vTfTsFv 
RveSd-l  ^ 

where  S=RTR  (u=Rv).  One  can  prove  ([4])  that,  if 
X-max(I  ®(-E+X.I)  +  (-E+X.I)®I  +  F®F)  <  0, 
then  the  equation 

S(-E+XI)  +  (-ET+X.I)S  +  F^SF  =  -I 

has  a  unique  solution  which  is  positive  definite.  Then,  for  any  X 
smaller  than  X.min(I  ®E  +  E®I  -  F®F)/2,  we  have 


lim  inf  |-log(IIXtll)  >  X 

t— >oo  Xq 

The  smallest  Lyapunov  number  of  Qg  is  then  larger  than 
^max(I  ®E+E®I-F®F)/2.  But 

I  ®E  +  E®I  -  F®F  =  -{  I  ®(A-B2)  +(A-B2)®I  +  B®B  }T 

This  ends  the  proof  of  the  theorem. 

We  have  just  proved  the  following  result 

Theorems :  The  smallest  Lyapunov  number  of  equation  (1)  is  larger 
than 


1 

6=2  ^niin(I  ®A  +  A® I  -  B®B). 

If  equation  (1)  is  replaced  by 

n 

dYs  =  AYs  ds  +  ^  BiYsdw^  (12) 

i=l 

where  Bi,...Bn,  are  n  matrices  and  wi,....,wa  are  n  independent 
Brownian  motions,  the  proof  can  still  be  carried  out  in  the  same  way 
and  we  obtain 

Theorem  1  (general  form):  If  equation  (12)  is  subsituted  to  equation 
(1),  theoreml  and  3  are  still  valid  with 
1  a 

Y=  r  ^max(I®A+A®I  I  ®Bi2+Bi2®I-Bi®Bi) 

^  i=l 

1  a 

5  =  “ A+ A® I  Bi®Bi) 

^  i=l 
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